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SEIiiPR FOR RADAR ASTRONOMY 

STANFORD UNIVERSITY 

DEPARTMENT OF ELECTRICAL ENGINEERING 

STANFORD. CALIFORNIA 94303 
September 13, 1983 


Henry C. Brinton 
Code EL-4 
NASA Headquarters 
Washington, DC 20546 

Dear Henry: 

The five attachments to this letter constitute the final report on 
contract NAGW-326, on "Continued Development of the Radio Science technique as 
a Tool for Planetary and Solar System Exploration". 

Attachment A introduces a possible alternative to a spacecraft monostatic 
radar system for surface studies of Titan. We believe that the potential of 
the uplink bistatic radar approach is so high that a detailed study should be 
conducted. We may propose such a study if there is a continuation of the PIDP 
program in 1984. 




Attachment B outlines the results of a short study of the characteristics 
of a bistatic radar investigation of Titan's surface, presented in terms of 
the Voyager 1 flyby and a proposed Galileo orbiter of Saturn. It shows the 
limitations o^ the downlink bi static approach and points the way to the need 
for an uplink experiment of greatly improved sensitivity. 


Attachment C is a general outline of the critical factors which need to 
be addressed in order to optimize che radio occultation technique for the 
study of clouds and cloud regions in planetary atmospheres. 

Attachment D addresses potential improvements in the techniques for 
measuring small-scale structures in planetary atmospheres and ionospheres 
The recent identification (Hinson and Tyler: 1983, Icarus . 54 , 337) of tne 

major role of internal gravity waves in Titan's atmosphere represents an 
important new departure in the general subject of small atmospheric structure, 
and “-he method used in isolating and identifying the gravity wave effects need 
to be perfected and applied to other atmospheres. 


Attachment E is a draft report on the development of a technique for # 
vastly improving the radial resolution from the radio occultation measurements 
of the rings of Saturn. While this work was supported primarily from other 
sources, the report is included here because extensions of this technique may 
well be applicable to obtaining increased height resolution in planetary 
atmospheres as measured by radio occultation, and improved areal resolution in 
planetary surface studies by the bistatic radar method. In all cases, the 
improvements would be due to, and would require, high levels of phase 
stability, as is demonstrated here from the use of the Voyager USD (ultra- 
stable oscillator) for the occultation studies of Saturn's rings. It is 
coiranonly said that the resolution is set by the size of the principal Fresnel 
zone, but in the examples here, resolutions to less than a hundreth of the 
Fresnel zone is achieved in the rings. 
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We believe that important progress has been made during the past year, 
based on the PIDP support. However, this is a continuing process. We suggest 
that follow-on efforts are needed in order to optimize the scientific 
potential of the radio systems that are continuing to be developed for 
spacecraft telemetry and tracking. 


Very truly yours. 



Von R. Eshleman 


63:Brinton3 
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STANFORD UNIVERSITY 

DEPARTMENT OF ELECTRICAL ENGINEERING 

STANFORD. CALIFORNIA 94305 
Septendser 7» 1983 


To: Titan-ites (CaTISiOg) 

From: Von R. Eshleman and G. Leonard Tyler, STAR LAB, 
Stanford University, Stanford, California 94305 

Re: Uplink Bi static Radar for Titan Surface Studies 


It is generally agreed that the next spacecraft mission to the Saturn 
system should attempt radar measurements of Titan's surface. Preliminary 
studies by European and U.S. investigators have included consideration 
of a spacecraft monostatic radar, such as the Pioneer- Venus Orbiter radar. 

A group at Stanford University is preparing a proposal to NASA for a 
detailed study of an alternate radar method, as outlined below. This memo 
is being sent to interested parties for their information relative to 
continuing planning of potential missions to the Saturn system. 

We suggest consideration of uplink bi static radar for detection and 
study of Titan's surface. While most bistatic radar and occultation 
investigations have been downlink (transmission from spacecraft and , 
reception on Earth), it is important here to realize the factor of 10'* to 
10^ improvement in sensitivity that should be realized with an uplink 
experiment using ground-based transmissions and spacecraft receivers. The 
"science Instrument" would be based on the use of the facilities of the 
NASA-JPL Deep Space Network (in particular, the large ground-based antennas 
and powerful command transmitters), the spacecraft telemetry antenna and 
feeds, and part of the spacecraft receivers. Provision of a low noise 
input stage and a wi de-band amplifier chain would be required on the 
spacecraft, with associated equipment for sampling and storage of data and 
for merging into the standard data stream. For the surface studies, the 
spacecraft high-gain telemetry antenna would be pointed toward Titan for 
the carom signals, but the direct signals would also be received, for 
reference, through the antenna sidelobes (or through a low-gain antenna, if 
one were already supplied for the telemetry system). That is, the only 
equipment supplied specifically for the experiment would be for amplifying, 
storing, and merging the wi de-band bi static radar data obtained by using 
the uplink command system. 
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The Voyager radio-science experiirenters considered downlink bistatic 
radar for Titan surface studies » but did not attempt it because even 
simple detection of a bounce signal appeared to be unlikely, and because 
it would have complicated spacecraft operations at a very critical time. 
However, the proposed uplink experiment would have orders-of-magnitude 
more sensitivity, so that it might well be capable of obtaining many of 
the Titan objectives of a spacecraft monostatic radar experin«nt similar 
to that of the Pioneer Venus Orbiter. Assuming a Saturn orbiter with 
repeated Titan encounters, as in the European Cassini study, the bistatic 
radar studies could include measurements of surface elevations (topography) 
along different tracks for each encounter, reflection coefficient and 
Brewster angle measurements to determine the aielectric constant of the 
surface material, spectral studies to investigate small-scale surface 
roughness along the radar tracks, and some degree of surface imaging 
based on tirr^ delay and spectral measurements. 


While it is clear that a monostatic spacecraft radar system could 
provide very significant data on Titan's surface, we believe that the 
uplink bistatic possibility should also be considered and investigated 
in detail, in case later study of trade-offs might make it a contender. 


If the uplink approach were chosen, the same system could be used 
for detailed investigation of Saturn's rings by the occultation 
technique, and repeated occultation measurements of the Titan and Saturn 
atirospheres , at even greater sensitivity and resolution than was obtained 
in the downlink Voyager expe*'i merits. That is, there would be continuing 
science to be conducted even if Titan's surface turns out to be 
everywhere the same; i.e., oceanic. 



Radar Scattering from Titan's Surface 
Using a Saturn Galileo Spacecraft 

Properties of a signal scattered from the surface of Titan during the 
Voyager I encounter were studied previously. Although some marginal 
detections appeared possible, no experiment was conducted. 

For Voyager, the most likely detection wwjld have occurred when the 
specular point velocity matched the body rotation - a sort of resonant 
condition. Even so, the signal had to be observed In a sub-Hertz bandwidth 
for several minutes before one could expect a detection. 

Similar calculations have been performed for a possible Saturn Galileo 
encounter with Titan. In fact, two cases have actually been considered - one 
using a spacecraft on a "polar tour" of the Saturnian system and one using an 
equatorial tour. Results for these, as well as the early Voyager study, are 
attached. 

The polar tour shows the least promise. Though closest approach is after 
occultatlon and total echo strength is highest then, so Is the bandwidth. 
Detectability - expressed as strength per unit bandwidth - is never larger 
than -45dB. 

For the Galileo equatorial tour there Is a grazing (perhaps atmospheric) 
occultatlon about 25 minutes after closest approach. Again, however, maximum 
strength and maximum bandwidth occur simultaneously. The bandwidth Is never 
less than 20 Hz for this trajectory so there is no possibility of a 
"resonance" detection. The best detectability for the equatorial tour Is 
about -44dB, shortly after closest approach but before deepest atmospheric 
penetration. 

The radio wavelength assumed for these calculations was 13 cm. 

Presumably use of 3,6 cm would offer higher transmit /receive antenna gains but 
with some additional loss to atmospheric absorption and (probably) a higher 
diffuse-to-specular scatter ratio. Highest detectabilities were obtained for 
grazing incidences, so these factors would have to be considered more 
carefully before any decisions were made. 

The surface of Titan was assumed to be perfectly reflecting with 6®rms 
(gently undulating) roughness. Lower reflectivity - which would be expected 
for a dielectric - would degrade the detectability. Increased roughness - 
such as is found on the outer three Galilean satellites - would also reduce 
chances for detection. On the other hand. Mars - another body with an 
atmospnere - appears to have only 2"rms roughness (and sometimes less) on its 
north polar cap. 

The "resonance" detection is very sensitive to Titan's pole vector and 
rotation rate - which, for this study were assimed to be perpendicular to Its 
orbital plane and synchronous, respectively. The Galileo trajectory is 
sufficiently oblique to Titan's equatorial plane that no resonance condition 
was found. 


To get the most out of a Titan BSR encounter, two conditions must be 
changed: 


a I The trajectory must be closer to Titan equatorial, and 
b) The encounter velocity must be reduced 

Given a detectable signal the following (partially dependent) paran^ters 
could be estimated: 

1) surface roughness 

2) surface dielectric constant 

3) pole orientation, and 

4) rotation rate 


Dick Simpson 
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PROGRESS IN DATA ANALYSIS; 

NEW TECHNIQUES FOR INTERPRETING RADIO SCINTILLATIONS 
DURING PLANETARY OCCULTATIONS 


D. P* Hinson, G. L. Tyler, and V, R. Eshleman 

The exploration of the Jovian and Saturnian systems using Voyagers 1 and 
2 included 3.6 cm and 13 cm wavelength radio observations during occultatlons 
of the spacecraft by Jupiter, Saturn, and Titan. One of several areas of 
research based on these observations Involves the way that refractive 
Irregularities In the atmospheres of these planets caused an^lltude and phase 
scintillations at both wavelengths. In conjunction with the reduction and 
analysis of these data (which was funded under a separate grant), we have 
developed new techniques for interpreting the scintillations observed during 
planetary occultatlons; each is described briefly below. In particular, 
there has been a significant iiTf)rovement in the procedure used to analyze 
scintillation spectra. Moreover, application of this methodology to the Titan 
occuUation data has yielded unique results which should help to resolve a 
long-standing controversy as to the nature of the atmospheric irregularities 
which cause scintillations during planetary occultatlons. 

Analysis of Scintillation Spectra. In this new procedure, measured 
temporal spectra of the weak fluctuations in amplitude and phase are 
interpreted using conventional methods of spectral data analysis and the well- 
developed theory for weak scattering. The occulting atmosphere Is modeled as 
a thin screen containing a power-law distribution of anisotropic 
irregularities. Comparisons of measured scintillation spectra with the theory 
produce least-squares solutions for four parameters which characterize the 


underlying atmospheric scattering structures. The parameters are: (1) the 

axial ratio* (2) the orientation, (3) the poMer-law exponent of the spatial 
spectrum, and (4) the magnitude of the spatial variations in refractivity. In 
order to determine the measurement uncertainties, the covariance matrix of the 
solutions for the four parameters is calculated by standard propagation of 
errors. 

Previous methods for analyzing the scintillation spectra relied on a 
simple visual comparison between the theory and the measurements, and provided 
no uncertainties in the results. The in^rovement offered by the new procedure 
is best illustrated through an example. We have analyzed the radio 
scintillations observed during the occultations of the Voyager spacecraft by 
the ionospheres of Jupiter (Hinson and Tyler, 1982, JGR, 87, 5275-5289) and 
Saturn (Hinson, 1983, JGR . in press). In particular, the magnetic field 
orientation at locations near these planets has been inferred from the 
occultation data, along with the associated uncertainties. These measurements 
are of interest because they provide a test of current Jovian and Saturnian 
magnetic field models, which were derived from spacecraft magnetometer 
measurements. In this regard, the uncertainties in the occultation 
measurements are of fundamental importance in evaluating the field models -- 
without uncertainties, the results derived from previous occultation studies 
are difficult to interpret. Hence, this new procedure has extended the 
capability of the radio occultation method. 

The Signature of Atmospheric Internal Gravity Waves in Occultation 
Measurements. In a wind-free isothermal atmosphere, the amplitude of internal 
gravity waves propagating without attenuation increases exponentially with 
altitude. F.eal planetary atmospheres have nonuniform temperature profiles; 
this thermal structure imposes additional variations in wave amplitude with 


height. The vertical wave envelopes can be calculated for each of the 
planetary atmospheres studies during spacecraft flybys on the basis of radio 
occultatlon measurements of tenperature-pressure profiles and a linearized 
theory for atmospheric waves. 

The density disturbances associated with atmospheric waves perturb the 
occultatlon radio signals. Under conditions of weak scattering, the 
fractional variance of received signal intensity is proportional to the 
variance in density. Hence, the presence of Internal gravity waves can be 
established through a convarlson of the theoretical calculations described 
above with the measured profile of the residual intensity variance (after the 
removal of average effects). In this manner, the signature of freely 
propagating internal gravity waves in Titan's upper atmosphere has been 
Identified in the Voyager 1 radio occultatlon data (Hinson and Tyler, 1983, 
Icarus, 54. 337-352). Above 25 km, the altitude dependence of the measured 
intensity variance appears to be explained entirely by modulation of internal 
gravity waves by the atmosphere's average vertical structure. 

Signal scintillations are a common feature of the observations during 
both radio and stellar occultations; this new method of analysis should lead 
to a better understanding of atmospheric effects in occultatlon 
measurements. In principle, this same procedure can be applied to the study 
of Venus, Mars, Jupiter, Saturn, Uranus, and Neptune. 


Clouds In planetary atmospheres 

I 

We have considered Idealized experimental approaches for using the radio 
occulatatlon technique to study cloud systems In planetary atmospheres. A 
summary account Is given here of our general conclusions concerning experiment 
design and types of measurements. As one result of this Investigation, we 
have prepared and submitted to NASA a proposal to use the Ploneer-Venus 
orblter for studies of the cloud region In the middle atmosphere of Venus 
(Stanford proposal number RL 38-83, submitted August, 1984, for Radio Science 
Experiments during the Pioneer Venus Extended Mission). That Is, while 
Pioneer Venus does not provide the Idealized radio science "Instrument" In the 
sense given above. It appears that It has significant potentialities of a type 
that have not yet been realized for this purpose. 

The general requirements for optimum radio propagation studies of the 
cloud regions In planetary atmospheres include: 

1. Highest feasible SNR (signal-to-nolse-ratio). 

2. Best possible frequency stability. 

3. Best possible amplitude stability. 

4. At least two and preferably more radio wavelengths. Including one sub- 
centimetrlc wavelength. 

5. Orthogonal signal polarization measurements including a pair at a sub- 
centimetric wavelength. 

6. Spacecraft orientation control for limb tracking and bistatic sounding 
of the surface. 

7. Dual spacecraft. 

Number (1) involves spacecraft and ground based transmitter powers, 
antenna sizes, receiver noise figures and bandwidths, and operating radio 




frequencies. There Is a long history of continuing efforts for Improvements 
In these areas by the Deep Space Net, In Its efforts to advance the state of 
the communications art for the radio links between distant spacecraft and the 
Earth. Essentially all in^rovements In these areas directly benefit radio 
science experiments. Item (2) has also been of prime concern to the DSN, but 
with emphasis only on the ground terminal. A new level of stability on the 
spacecraft was achieved with the USD (the so-called ultra stable oscillator) 
Introduced for radio science on the Voyager missions. The USD constitutes an 
extraordinary exanqile of the value of a modification to the communications and 
tracking systems for the purpose of Improving radio science. Fundamentally 
Important aspects of the atmospheric and ring occultatlon experiments, and 
certain tests of the theory of relativity, were possible with Voyager only 
because of the USD. Still greater spacecraft oscillator stability Is possible 
for future missions, and would directly benefit measurements of atmospheric 
structure at precisions that would be Important In helping to determine the 
effects of condensibles. 

Amplitude stability (Item 3} has not been an area that has been 
enphasized (except for the usual high standards of engineering followed by the 
DSN), since there are no stringent requirements for the communications and 
tracking applications. For the study of cloud regions and their associated 
absorbing vapors, htwever, amplitude stability factors directly Into the 
accuracy of the absorptivity measurements. Thus this Is a key requirement for 
cloud studies. The ability to measure small amplitude changes with accuracy 
Involves not only the knowledge of transmitter powers and receiver gains, but 
also antenna gains In the propagation directions for both the spacecraft and 
ground-based antennas. This Implies stringent requirements on the knowledge 
of the pointing of these high gain antennas. There Is much to be gained In 


the potential of future cloud studies from a concerted effort to Improve our 
knowledge and control of signal amplitudes. 

The history of the debate over the need for multiple radio frequencies or 
wavelengths In spacecraft radio experiments (Item 4) goes back to the time of 
the first deep-space missions. The Joint advantages of the S and X band 
systems for communications^ tracking, and science resulted in their eventual 
development, and has thus far led to their continued use. For cloud studies, 
multiple wavelengths are Important on at least four counts. First, they allow 
self-calibration of plasma effects on the long links to Earth, for greater 
profile accuracies. Second, the dependence of absortivlty on wavelength can 
help Identify and characterize the cloud vapors and possibly the particulate 
absorbers. Third, signal scintillations (amplitude and phase) at several 
wavelengths give a more con^lete picture of small-scale atmospheric structures 
that may be related to the effects of the latent heat on the cloud 
condensibles, convective Instabilities, and Internal gravity waves. And 
fourth, the use of a relatively short wavelength may make It possible to 
measure the depolarizing properties of large cloud particulates, as discussed 
next. 

Depolarization has been measured on radio propagation paths between 
communications satellites and ground tracking stations. Moreover, the 
character of the depolarization has been observed, on occasion, to change 
suddenly. It 1s believed that the depolarization Is due to Ice crystals and 
that the changes accompany disorientation or reorientation of crystals that 
Initially were aligned by atmospheric elactric fields. That Is, the sudden 
changes have been observed to accompany lightning discharges. Thus It Is 
possible that polarization measurements on occultation paths through the 
atmospheres of the planets (item 5) would provide information on the number 


and size of certain cloud particulates, about the^r shapes and orientations, 
and about sudden changes In electric fields In these atmospheres. The 
emphasis on relatively short wavelengths for such measurements Is because such 
effects would then be accentuated. 

Spacecraft orientation control (Item 6) Is related to the pointing 
accuracy problem as described for Item 3 on amplitude stability. In addition. 
It may be possible to obtain limited data on cloud related properties for 
regions below the level of critical refraction by a bistatic radar bounce off 
the surface. This Is particularly relevant to Venus, where the occultatlon 
measurements can probe down only to about 33 km altitude. Much of the region 
of Interest for Venus cloud studies lies above this height, but aborptlon by 
cloud related gases certainly occurs In the lov^r regions that can only be 
reached by such a radar bounce. For this experimental approach. It would be 
necessary to stu4y also the scattering properties of the surface, since they 
must be known to determine the atmospheric absorption along the carom path. 

Two or more spacecraft slimiltaneously In orbit around the target planet 
(Item 7} would provide. In principle, two methods for In^rovlng the accuracy 
of occultatlon experiments. Consider one spacecraft In atmospheric 
occultatlon and the other well away from the planet In the plane of the sky. 
Simultaneous radio links from both to the Earth would provide atmospheric 
occultatlon signals, plus "calibrating" signals that traverse essentially the 
same path In Interplanetary space and In Earth's atmosphere and Ionosphere. 
Thus the effects of the planetary atmosphere could be Isolated more thoroughly 
than In the one-spacecraft method. Tie other possibility for multiple 
spacecraft Is for the occultatlon radio link to be between pairs of orbiting 
spacecraft. The signal charasterl sties could then be measured with the 
Information sent to Earth by telemetry. Because these occultatlon paths would 




be very short compared to the distance to Earth, higher SNR and less 
contamination by terrestrial and Interplanetary media could be obtained. In 
principle. 

It would of course require considerable time, money, and effort to 
address all of the areas where significant Improvements could be made In radio 
occultatlon studies of planetary atmospheres, with the above en^ihasls on the 
regions of clouds. It took many years to progress from single to dual 
frequency communications and tracking systems, and the USO Is the only major 
example where a modification of these systems was made primarily for radio 
science. Thus It seems Important to Identify an order of priority. Assuming 
the continuance of dual frequency systems (4), we believe that the next order 
of priority Is Increased attention to Improving the accuracy and reliability 
of amplitude measurements (3 and 6), and further in^rovement In frequency 
stability (2) for the spacecraft source beyond that presently provided by the 
USO or equivalent. The introduction of at least a third radio wavelength 
(shorter for some purposes and longer for others) would be next on our 
suggested agenda (4), with other proposed In^rovements being Implemented when 
feasible. 


FINE-RESOLUTION RECONSTRUCTION OF THE MICROWAVE OPACITY 
OF SATURN'S RINGS: THEORY. EXAffLES, AND LIMITATIONS 


Essam A. Marouf 
G, Leonard Tyler 
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Abstract 


INTRODUCTION 

Diffraction of radio waves Is a prominent phencMnena In measurements 
obtained during Voyager 1 radio occultatlon by Saturn's rings (Narouf and 
Tyler, It Is caused by variations In the microwave opacity of the 

rings (optical depth at radio wavelength) on a radial scale comparable to or 
smaller than a characteristic dimension whose value is determined by the 
occultatlon geometry and observation wavelength (Fresnel zone radius F). For 
Voyager 1 geometry, 9 £ F 15 km at 3.6 wavelength (X.) and Is nearly twice 
as large at 13 cm-X, setting a resolution limit of few tens of kilometers on 
opacity profiles obtained without corrections for diffraction. Although such 
resolution Is certainly adequate for revealing a great deal of structure 
within the rings (Tyler et al. , 1983), a much finer resolution becomes 
essential If the opacity profiles of the highly localucd ring features 
discovered by Voyager (e.g. , sharp edges, narrow ringlets and gaps, wave 
phenomenon) are to be properly mapped. A major objective of this paper Is to 
develop Inversion algorithms that correct for diffraction effects, hence 
Improve achlevalbe resolutions to values « F. For Voyager 1 at Saturn, tne 
refined resolution Is smaller than F by a factor « 10 to 100, depending on 
the opacity of the region probed and the measurement errors. If Voyager 2 
completes Its planned occultatlon by Uranus* narrow rings successfully 
(January, 1986), a Fresnel zone radius ■ 2 kin at 3.6 cm-X will make the need 
for diffraction corrections almost Indispensable. 

In radio occultatlon of planetary rings, opacity Information Is obtained 
from measurements of the average effects cf ring material on an Incident 
monochr(Mnat1c plane radio wave traveling along the 1 1ne-of-s1ght between a 


2 


spacecraft and the Earth. By utilizing an ultrastable oscillator to generate 
this wave, (Eshleman et a1, . 1977), not only the extinction of the direct wave 
emerging from the rings can be measured, but also its phase relative to that 
of the Incident wav-^ (hence the terminology "coherent signal"; Marouf et al. . 
1982). It Is the abMIty to simultaneously measure the amplitude and phase of 
the coherent signal that makes the fine resolution reconstruction of the 
opacity profile possible. Such measurement correspond. In fact, to a one* 
dimensional sampling of the microwave hologram formed behind the rings. 
Starting from Fresnel -Huygens formulation of the diffraction problem and using 
a planar axially symiwtric screen model for the rings, a Fresnel transform 
relationship between the ccm^ex amplitude measured (the hologram) and the 
microwave transmittance of the ring-model Is shown to exist. The 
corresponding Inverse transform provides the desired transmittance (hence, 
opacity) with resolution « F. Factors analyzed that affect post-inversion 
resolution are: uncertainties In F, finite Cv.'rvature of ring contours, 

reference osclllaotr stability, and random measurement noise. Although the 
methods developed are general enough to be appllcaf^e at Uranus and also to be 
useful In guiding choices of optimum occultatlon geometries of future 
missions, our emphasis throughout this paper is on their application to the 
Voyager 1 occultatlon by Saturn's rings, where the geometry and selecte^d 
examples of the measurements are used to Illustrate the Inversion theory and 
its limitations. Other examples are given by Tyler et al. (1983). A more 
complete diffraction corrected opacity profiles over the full extent of theX 
rings at 3.6 and 13 cm-X Mill subsequently be available. 
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INVERSION THEORY 


Our objective In the analysis to follow Is to show that there exists a 
Fresnel transform relationship between the (complex) amplitude of the observed 
coherent signal and the microwave transmittance of the rings (defined 
below). Consequently, single analytic Inversion of such transform allows 
reconstruction of the transmittance at a much finer resolution than would have 
been possible had the phase of the coherent signal not been also measured. 
Geometry 

The geometry of radio occultation and formulation of the received signal 
in terms of physical ring models have been discussed at length elsewhere 
(Marouf et al. 1982) and results of diffraction limited microwave opacity 
profiles over the full radial extent of the rings have also been published 
(Tyler et al. , 1983). Fig. 1 gives an alternative view of the occultation 
geometry projected on two orthogonal planes in a Saturniocentric coordinate 
system. The xy-plane coincides with the mean ring plane, which is assumed to 
be orthogonal to the north pole of Saturn (z axis). During occultation, the 
spacecraft and Earth are located on opposite sides of the rings, hence the 
line-of-sight penetrates the ring plane at the point (Xq, y^), or (p^, 4 »^) 
if circular cylindrical coordinates are used. Azimuthal angles are referenced 
to the x-axis whose direction is along the Earth-Saturn vector at the time of 
the experiment projected on the ring plane. In Fig. 1, the angle ^defines 

A 

the ring opening, is a unit vector along the line-of-sight, and R^. is 
the spacecraft position vector; other geometric parameters are as defined by 
the Figure. Typical values of some geometric parameters for the Voyager 1 
radio occupation experiment (November 13, 1980) are given in Table I; other 
entries to this Table will be introduced later. 


4 


Please note that in the actual experiment, dual-wavelength coherent 
signals were transmitted from the spacecraft and received on the Earth. 

Because it is somewhat easier to think in terms of an incident plane wave and 
diffracted spherical waves, the analysis here assumes the opposite 
configuration of the transmitter and the receiver; by reciptocity, however, 
identical results should be obtained for the actual experimental 
configuration. 

Fresnel -Huygens Formulation 

The reader is reminded here that because of the random positions of ring 
particles in space, and also because of the random intrinsic properties of the 
particles (e.g. , size, shape, orientation, ...etc.), the radio signal emerging 
from the rings during occultation is random in nature. This paper deals only 
with the coherent component of this signal, that is, the con^onent that 
characterizes the average effect of ring material on the ray directly 
transmitted through the rings (Marouf et al. , 1982). Such effect can best be 
described in terms of the Fresnel-Huygens principle as folltws (v. , e.g., 
Goodman, 1968). 

One can think of the average effect of the ring material on the incident 
wave in terms of the effect of an extended planar "screen" which coincides 
with xy-plane in Fig. 1. Thus, for a plane harmonic wave of frequency w, 

A 

wavelength X, and waver.umber k (k=2u/X) indIdent along u^, the coherent 
field observed at R^. may be obtained by coherent addition of the 
contributions of all (fictitious) secondary sources formed just behind this 
screen. The result is 


iku.'p 

= J/ dp T(p) e ’ ] 


e'ElVI 


( 1 ) 
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where the term In brackets represents the complex amplitude of a secondary 
source of infinitesimal area dp and the last term repres(*nts the spherical 
wave that eminates from such source to the receiver at Note that the 
complex amplitude Is affected by: (1) the so-called "obliquity factor," which, 
because the occultatlon geonwtry limits observation angles to near forward 
directions. Is asswned in Eq. (1) to be constant ■ « cos(B), (11) the 

relative phase of the Incident plane wave at p, represented by the term 

A 

exp(1ku^*p), and (111) the deterministic function T(p) which embodies the 
average amplitude and phase changes Introduced by the ring particles. One can 
define T(p), therefore, as the ratio of the complex amplitude of the 
secondary source formed behind the ring plane at position p In the presence 
and absence of the rings, respectively. We refer to T(p) in this paper as 
the (complex) transmittance of the rings. Our objective will be to use the 
measurements of E^. to recover T(p) from Eq. (1). 

Before we proceed to acconpllsh such objective, we introduce the new 

A 

normalized observable function T defined as 

A 

A -iku.'R 

T - E^e ’ ' (2) 

that Is, the observed coherent field normalized by the incident field that 
would have been observed at R^. had there been no ring material along the 
1 ine-of-sight. Such normalization compensates the phase of E^. for a geometric 
component introduced solely by the changing phase path due to spacecraft 
motion along its accurately known trajectory (Tyler et al. , 1983). Thus both 

A 

T and T embody only effects caused by encountered ring material. 

The question now is: can one reconstruct from the measurements of the 

A 

complex samples of T versus time the transmittance function T(p) from Eq. 


6 


(1)? The answer depends, of course, on wrethar further assianptlons can be 
made about T(p). In general, the two dimentlonsi function T(p) cannot be 

A 

reconstructed from measurements of T a'.ong a single occultatlon track. 
Fortunately, however, the near circular symmetry of Saturn's rings, at least 
over the limited region of ring area effectively contributing to Eq. (1), 
renders T(p) to be one dimensional (function of the magnitude p only) thus 
allowing for a positive answer to the above question. 


ID-Screen Approximation 

In a circular cylindrical t.oordinate system in which the screen coincides 
with the (p,^) plane and T(:; = T(p), Eqs. (1) and (2) combine to give 


T 



2 % 

dp p T(p) / d« 
0 


gik[|R^- p|-(R^-p)*u^] 



( 3 ) 


Note that when p = p , R -p* is In the direction of (See Fig. 1) and the 
exponent in the ♦-integral vanishes. At positions p far away from this 
(stationary) position, the ♦ integrand oscillates rapidly, resulting in 

A 

negligible contribution to T. The extent of the screen area that 

A 

significantly contributes to T can be found by assuming ♦ = ♦, where 

♦ is small, and performoing a perturbation expansion in powers of ♦. 

Details are given in the Appendix; the result when keeping only first order 

Mr 

terms in ♦ is 


T(Pq) 


1-il 



dp/F 


( 4 ) 


where the characteristic dimension F (the interval P - Pq for which the 
Kernel phase = n/2) defines the radius of an effective Fresnel zone along the 
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radial dimension p, whose value Is solely determined by the occultatlon 
geometry and observation wavelength X, namely. 


F 


2 2 2 
M) S ^0 * ^^jl/2 


T 


( 5 ) 


where D » *^o “ * cos ♦q. » sin ♦q, and the angles 

B and have been defined in Fig, 1, Typical values of F computed for the 

Voyager 1 radio occultatlon geometry at x ® 3,6 cm are given In Table I, 

The Infinite limits of the Integral In Eq, (4) are, of course, only a 
mathematical convenience, since the actual extent of the ring material Is 
finite In p. These limits allow one, however, to Interpret Eq, (4) as a 

A 

Fresnel transform relationship between T and T, whose Inverse can be shown to 
be 


T(P) = 1 T(p„) e 


'"'’ 0,2 

F 


dPo/F 


( 6 ) 


Thus, from the occultatlon geometry, the Fresnel zone F Is first computed, 

A 

then the measured complex samples of T are "filtered" using Eq. (6) to 
reconstruct T at radial resolution much finer than F Itself, An example 
to demonstrate the working of the theory and to motivate the discussion to 
follow of Its limitations Is in order. 


Example I 

Numerical computations of the inverse Fresnel transform (6) can be 
carried out using either one of two approaches. The first is a 
straightforward discretization of the integral using the measured temporal 

A 

samples of T interpolated on an equispaced grid in p^, namely. 
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T 


m 



N/2 . 


-I7 Tl*- 

/ 1 e dTi 

(nt^)Ap/F 


( 7 ) 


where Ap is the distance betv^en samples of T, is T(p » mAp), m an 
integer, and NAp is the width of the data window (centered at mAp) used to 
evaluate The Fresnel integrals on the right side of (7) can be 

computed, for example, fr(xn known series approximations (v. , e.g. , Boersma, 

J., 1960). 

The second approach relies on recognizing the fact that the integral in 
Eq. (6) is actually a convolution integral which may be evaluated using Fast 
Fourier Transform (FFT) techniques. Because of its speed, the FFT approach is 
more convenient for systematic reduction of the large voltmie of data collected 
over the full extent of the ring system. Both approaches were used in our 
data reduction; with excellent agreement obtained. 

Experimental details relevent to the radio occultation data reduction 
procedure have been discussed at length by Tyler et al. (1983) and will not be 
repeated here. We draw the attention of the reader, however, to the fact that 

A 

T in the analysis here is equivalent to E in the notation of Tyler et al. 
normalized to unit amplitude in free-space. 

Our first example is derived from the Voyager 1 data in the neighborhood 
of the Encke division, whose center is located at p* 133,590 km and whose 
width » 320 km. Both edges of this division exhibit clear diffraction patterns 
that are characteristic of two widely separated knife edges of gray screens 
(Marouf and Tyler, 1982). Fig. 2a and b show the measured magnitue and phase 

A 

of T generated from independent data samples separated a distance 

Ap = 1.2 km. Compared to a Fresnel zone radius of F = 15,273 km at the mean 
radial position of the division, this separation distance adequately samples 
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the first several oscillations on the clear sides of the edges. To obtain the 
inverse transform, Eq. (6) was used, assiming F above and a data window of 
width W * NAp = 390 m. The mcgnitude and phase of the recovered transmittance 
T{p) are shown in Fig, 2c and d, respectively. Note in particular the 
disappearance of any regular oscillation pattern in the magnitude curve, the 
absence of any marked overshoot of free space level, and the much sharper 
signal level transition at the edges. Numerical examination of sample values 
of the reconstructed magnitude indicate that the transition from near-zero to 
near-unity levels occurs over the distance of one sample at either edge, in 
agreement with an upper bound of about 1,5 km on edge taper projected in the 
ring plane obtained by model fitting studies of the same observations (Marouf 
and Tyler, 1982), 

A first order conclusion based on the results above is that the inversion 
algorithm appears to work. However, several important issues remain to be 
examined before one can draw any quantitative conclusions from the recovered 
profiles. Obvious questions are: (i) does the data window width W chosen 
preserve the resolution of the original data points?; (ii) can the resolution 
be further refined so that tighter bounds on the edge taper may be obtained?; 
(iii) what are the limits on the resolutions achievable in other thicker parts 
of rings?; (iv) how sensitive is the inversion process to perturbations in F?; 
(v) does any of the signal fluctuations within the apparently clear division 
(see Fig. 2c) correspond to a real embedded feature?; and (vi) when does the 
ID-screen approximation used become invalid? 

To shed more light on these important issues it is necessary to examine 
in some detail the systematic and random limitations of the Inversion theory 
above. This is carried out throughout the remainder of this paper, where the 
above example is revisited and other examples are also considered. 
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SYSTEMATIC LIMITATIONS 


Much like Fourier transforms, the forward and inverse Fresnel transforms 
(Eqs. (4) and (6)) are symmetric in form, except for the sign of i. One may 
think of the forward transform as an operation that distorts the actual 
profile T(p) through convolution with quadratic phase factor, and of the 
inverse transform as the “filter" necessary to compensate exactly for such 
distrition, so that T(p) may be recovered. If the inversion filter is 
slightly imperfect, either because a finite data window W is used and/or 
because F is not exactly known or vary across W, the reconstructed profile 
will be a somewhat distorted version of T(p). Even in the absence of any 
filter in^erfections, errors in T(p) can be caused by systematic drifts in 
the phase of the reference oscillator used. Below, we quantify the effect of 
these imperfections in terms of a single parameters, namely, the radial 
resolution R, and examine conditions under which the simple inverse 
transform (6) needs extension. 

In each of the cases examined, we denote by H^yg(p) the distorted 
profile recovered with an imperfect filter from the "data" corresponding to a 
hypothetical transmittance T{p) = 6(p), where 6(») is the Oirac-delta 
function and the subscript SUB identifies the source of the imperfection. One 
measure of the radial resolution Rjyg is the finite width of Hgyg, taken 
here to be the equivalent width of iHgygl^ (v. , e.g. , Bracewell, 1977), 
namely. 




'SUB 


SUB' 


( 8 ) 
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Our objective will be to find which of the imperfections considered sets the 
limit on the finest resolution achievable and determine such limits. 


Finite Data Window 

In principle, the diffraction pattern of any local ring feature extends 
indefinitely in space. Thus by limiting the extent of the data to be inverted 
to a finite width W, one could only hope to recover the transmittance of 
such feature to within some finite resolution. Think, for example, in terms 
of the classical diffraction pattern of an opaque knife edge. On the clear 
side of the edge, assumed to be at p^, the period of oscillation of the 
magnitude of the diffraction pattern decreases with distance from the edge 
approximately as 2F /|p - p^j. By chopping off the pattern everywhere 
outside the interval |p - p^[ i. W/2, all spatial frequencies of wavelength 
i. 4 f2/w are lost, thus when this "data" is inverted, the edge must exhibit a 
finite rise distance whose magnitude is of order 4F*-/W. 

As outlined above, we approach the resolution problem quantitatively by 
computing first the distorted profile H^^(p) corresponding to T(p) = 6{p). 

A 

Substituting this impulsive profile in Eq. (4), using the resultant T(p^) 
back in Eq. (6), and performing the integration over the finite Interval 
-W/2 1 Pq i. W/2, one gets 


Hw(p) 


1-^ s1nct"^ p)i e 

2F^ 2F"^ 


i ^4 
2 F^ 


( 9 ) 


where sinc(x) = s1n(xx)/iix. The original impulse T(p) and the magnitude 
lHy(p)l shown in Fig. 3a and b, respectively, plotted as a function of 

the normalized distance p/F for the case W/F * 50. As one may expect, a 
finite W maps the original 6-function into a new profile which possesses a 
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main lobe of finite width and many other minor side lobes. The equivalent 
width of sucn a profile Is obtained by substituting £q. (9) In Eq. (8); the 
result Is 


Ry - 2P2/W (10) 

a relationship that Implies that the resolution can be made arbitrarily fine 
by simply increasing W. This 1s true only If the Fresnel zone F Is exactly 
known and remains constant over W, an obvious Idealization that Is removed 
below. 


Fresnel Zone Uncertainty 

Substantial distortion In the shape of H^(p) "'ay result if the Fresnel 
zone F used In the Inverse transform (6) Is In error by an amount 6F. 
Generalization of Eq. (10) to include the effect of 6F Is carried out In the 
Appendix, where we show that If 16F/F| « 1 then Eq. (10) assimies the 
modified form 


H^p(p) = CC(5j) - C(^) + iS(5i) - iS(52)3 


-i ^ (■ 


7 


W 


(11a) 


where y = /26F/F and C(») and $(•) denote the real and imaginary parts 
of the Fresnel integral (Abramowitz and Stegun, 1972, p, 300). The arguments 
and ^2 are given by 


' yT 


7F 


(lib) 
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2 

An example that demonstrates the marked sensitivity of |H^p(p)| to 
errors In F Is shown In Fig, 3c, where 6F/F ■ 0.3% and W/F ■ 50, the same 
window width used In rig. 2b where 6F « 0. Note In particular the 
substantial distortion in the shape of the profile and the marked Increase In 
Its width for an error as small as 0.3% In F. 

The radial resolution corresponding to can be obtained from 

Eqs. (8) and (11) ^nd 1s shown In the Appendix to be 


R 


6F 


2F' 


W C^(Q) + S^Q) 


(12) 


where Q ■ yW/ 2F. In the limit tW 0, Eq. (12) reduces to Eq. (10), as It 
should. Contrary to R^ In Eq. (10), however, ncreases linearly 

with W in the limit W ♦ •. Salient characteristics of R^p are displayed 
by the solid curve in Fig, 4, where R^p/yF is plotted versus Q. The dashed 
lines represent the asymptotes yW 0 and yW •. The optimal resolution 
Is achieved when Q “ 1, that Is, when the width of the data window Is 


“6F ■ 


(13) 


and has a value 


- 2.5 F^/W^ » 1.25 yF (14) 

which is 25% larger than the 6F * 0 value obtained from (10) for the same 
window width. 

Eq. (14) implies that to achieve a sub-kilometer resolution 

3 

(RgP < 10 m) In regions where F Is 9-15 km (Table I) requires y to be 
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i 0.05, or |^/F| 4 0.14%; sub-hundred meter resolution, on the other 
hand, requires that |dF/F| ^0.0014% 

This marked sensitivity of resolution to errors In F Irameulately 
prompts the question: Hom accurate can F be estimated from the gaometry of 
the radio occultatlon experiment? As one may see from Eq. (5) for F, the 
answer de( jnds on the accuracy In knowledge of the geometric parameters D, 
and Mq* SucN parameters are affected by uncertainties In knowledge of 
the direction of Saturn's pole vector and also by timing uncertainties In the 
position vector of the spacecraft along Its trajectory. These uncertainties 
primarily translate Into an error 6p^ In determination of the radius of 
Intersection (Fig. 1), from which 6F can be approximately estimated to 
be GF ■F'6p^, where F' ■ bF/dp^ evaluated at the nominal point of 
Intersection. 

Utilizing the different viewing geometries of the radio and ultraviolet 
Voyager occultatlon experiments to correct for the offset In the radial 
locations of several prominant ring features, Simpson et al. (1983) report 
corrections to the direction of Saturn's pole that constrain 6p^ to be well 
within 10 km. Conservative bounds on i6F/F| in different regions of the 
rings can thus be obtained from the slopes F' listed In Table I using 
|6 pq| i. 10 km; such bounds were used to compute the values of y given In 
Table I. Corresponding optimal resolutions and window widths can be computed 
from Eqs. (r.j and (13), respectively, and are given in Table II, where *^6F 
varies from a few hundred meters in Inner ring C to a few tens of meters In 
the neighborhood of ring F, Whether these resolutions are achievable cr not 
depends on the limitations imposed by tiie additional Imperfections discussed 
below. 
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Finite Curvature 

One Idealization In deriving the Inverse transform (6) Is the asstanptlon 
of a constant Fresnel zone F. Two effects combine to violate this 
assiaaptlon: The first Is the time changing occultatlon geometry which 
causes F to vary with p^. as the finite, yet small, derivative F' In 
Table I Indicates; the second Is not as obvious, and Is caused by the finite 
curvature of the assumed circular contours of T(p). The analysis below shows 
the second effect to be d«n1nant for the Voyager geometry. 

An alternative approach to derive the forward transform (4) Is to 
approximate T(p) at any observation point by a new function which 
assumes constant values along straight lines tangent to t.he contours of T(p) 
and perpendicular to p^. As the width W Increases to accomodate ' 
resolutions, this straight-contours approximation becomes less valid ,.j 
( 4) must be extended to Incorporate at least first order effects or the finite 
curvature. 

Recall that Eq. (4) was obtained from the general form (3) through the 
assumption that over the ring area effectively contributing to the coherent 
signal the angle p deviates a small amount 4 from 4^, so that a 
perturbation expansion In 4 leads to Eq. (4) If only first order terms In 
4 are kept. This approach allows one to extend Eq. (4) to Include tno first 
order effects of curvature through Incorporation of terms upto second order In 
4 In the perturbation expansion. In the Appendix, we show that ti.e dominant 
effect of the additional terms In 4 Is to modify Eq. (4) to the new form 


T(Po) 


1-11 


J T(p) e 


1 


2F! 


-(P - Po)J 


dp/F 


(15) 
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where Is a dimensionless quantity that plays the same role as a partial 
derivative of F with respect to p evaluated at the point of observation 
Pjj. Its ma^ltude Is inversely proportional to the radius of curvature and 
Is given by 


F' 


' ■ tr— (u 


-1,2 ,2,2 
'U* ) c„s„ 
0 ^ 0 0 


+ s 


0*^0 g 

1 2'7'T 
s u + c 
0*^0 0 


(16) 


where the geometric parameters Cq» and Sq have been previously 

defined* Table I gives typical values of f; for the Voyager geonetry, 
cort^jted from Eq. (16) using fable I entries for F, and p^. Note 
that F^ as defined by Eq. (16) Is always positive. 

Contrary to the forward transform (4), the modified transforr (15) does 
not seem to possess an exact analytic Inverse similar to Eq. (6). One 1s 
ten^ted, however* to still adopt the syiranetrlc form 


P-P, 


2F- 


^(p) 


1+11 


/ T{p 1 
PlW/2 ® 


(P-P^)j 


dPo/F 


(17) 


as an approximate Inverse of Eq. (15). Clearly, T(p) will differ from the 
original profile T(p), and such difference will set an additional limit on 
achievable resolution. Before we compute uch limit. It Is Instructive to 
examine first the additional effect of the geometric change In F. 

Generalization of Eqs. (15) and (17) to Incorporate the geometric change 
In F Is heurlstlcally obtained by substituting the variable Fresnel zone 
F(Pq) for F. For the Voyager geometry, the drift In F with Is 
dominated by a linear component and can be locally approximated by the first 
order Taylor expansion 


17 


(18) 


f(Po) * F(P) + (P„-P)f 


where typical values of F* are given In Table I. If, In addition, the 

relative perturbation In F Is small (I.e., If ^)» 

first order perturbation In the phase of the Kernel of Eq. (17) is of 


Importance, giving 


P±W/2 ° ° 


as the modified approximate transform that compensates for both geometric and 
curvature-caused drafts in F, where F now actually stands for F(p). 
Interestingly, since F* and F^ are of the same sign (see Table I), their 
effects in Eq. (19) act to cancel one anthoer. Unfortunately, however, their 
magnitudes greatly differ, with the effect of F* being the dominant one. 
Thus, in the following discussion of resolution we limit our attention to Eqs. 
(15) and (17) only. 

As before, we obtain the resolution from the profile T(p) = U^(p) 

corresponding to T(p) = 6(p - p*), where the subscript C stands for 

curvature. Contrary to the assumption p* = 0 before, p* here is assumed 

large so as to avoid a fictitious singularity in F' since p assumes 

values of the same order as p"* (see Eq. (16)). The profile Hq is then 

obtained by substituting Eq. (15) back into Eq. (17) and performing the 

W F* 

resultant integration. The result (for 7 Tp— « 1) is exactly similar to 

I “ • I 

Eq. (11) except that the constant y in. that equation. Is nov< replaced hy the 
p-dependent term [6F^(p - p*)/F]^^^. The corresponding equivalent width is 
obtained from Eq. (8) and is shown in the Appendix to be 
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( 20 ) 


where U = 3F^ W/F and is assumed to be < 1. Thus, In the limit 0 

(I.e. , U > 0), Eq. (20) reduces to Eq. (10), while in the limit U 1, Rq 
increases indefinitely (probably an artifact of the violation of the initial 
assumption jF^W/2F| « 1), A minimum of Rq is reached at U » 0,8, as 
illustrated by the dash-dot curve in Fig. 4, where Rq/ 3F^F is plotted 
versus U. The corresponding optimal width Wq and resolution Rq are given 
by 


Wj. » 0.27 F/F; (21) 

and 

“ 2.75 F^/W^ = 10.30 F^F (22) 

A comparison of values of Rq and Wq to R^p and W^p is given in 
Table II. Clearly, Wq and W^p are comparable in ring C, both limiting the 
resolution to few hundred meters. The two widths systematically diverge, 
however, with increasing distance from Saturn, with Wq becoming the dominant 
limiting factor, where it limits resolution to about 100 meters in outer ring 
A and in the neighborhood of ring F. 

Oscillator Phase Drift 

Recall that the inverse transform (6), or its modified form (17), 
requires full knowledge of the amplitude and phase of the measured signal 

A 

T(p^). The phase measurement is actually a relative measurement, referenced 
to the phase of a coherent sinusoidal signal generated onboard the spacecraft 
from a highly stable, radiation hardened, quartz oscillator (Eshleman et al. , 
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1977), The radio occultation measurements were preprocessed so as to 
compensate for large scale systematic drifts in the estimated frequency of 
this oscillator (Tyler et al, , 1983). Nonetheless, a residual phase error 
♦j,(Po) nwy persist over the width of the data window used and the effect of 
these residuals on the inversion is assessed here. 

Consider the inverse transform (6) in the presence of slowly drifting 

A 

residual phase ♦p(Pq) contaminating the phase measurements of T(p^). 
Assuming that may adequately be approximated by the first three terms of 
its Taylor series expansion about the center p of the data window, i.e., 
V^^o^ “ ^p(p) ^r '^^0 ” perturbed inverse 

transform beconres 


T(P) = 


1+il 

~r~ 


i<tp(p)-i«p(p) 


i T(pjj)e 


HlP«-i 


r'^o 






dPo/F 


(23) 

It is easier to assess the effect of each of the expansion terms of 4 >p(Pq) 
separately. Clearly, the first term ♦^,(p) is a simple multiplicative 
constant. In the absence of the linear phase drift in (23) introduces 

a linear shift in T(p), whose value can be calculated by substituting in 

A 

(23) the profile T(p^) corresponding to T(p) = 6(p). The result is 

T(p) - 6tp + e (24) 

2 

Tnus, 4>p primarily introduces a shift = /n in the radial positions of 
the recovered profile. On the other hand, in the absence of <l>^, the 
term in Eq. (23) mimics the presence of an error 6F of magnitude 
6F/F = ()>^'F /27t, whose effect is to distort the recovered profile, as 
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discussed before. In general, a combination of shift and distortion would be 
expected. 

Estimates of and for Voyager radio occultatfon were obtained 
from polynomial least-mean-square fit to the residual phase measured in free- 
space before and after Voyager was occulted by the rings. Several data 
windows of typical width in the range specified by Table II were used, leading 
to the bound 10“^ rad/km, that is, a radial shift uncertainty i 70 

meters for F = 9 to 15 km. Estimates of fall in the range 

l^p'l ~ rad/km^, or an apparant error 6F/F £ , which has 

insignificant effect on values of y listed in Table I. One can conclude 
therefore that for Voyager conditions at Saturn residual oscillator drift have 
little effect on resolution and mainly introduce uncertainty in the radial 
position of value less than several tens of meters. 

Example I Revisited 

A comment about the resolution of the inverted profiles of Example I 
above is in order. For that example, F = 15.273 km at p = 133,590 km; 

other relevent parameters are: W = 390 km, F' = 1.6 x 10"^, 

Y = 0.0047, F' = 7.4 x 10"^. . Using these parameters in (13) and (21), 

one obtains W^p = 6500 km and Wq = 5573 km. The last value sets the limit 
on achievable resolution to Rq = 116 m. Because W << Wq, the resolution 
corresponding to W used is given by Eq. (10) to be = 1.2 km, which is 
actually the distance Ap^ used to sample the initial diffraction limited 
data in Figs. 2a and b. Clearly, ample room for further improvement in 
resolution remains. The example is considered again at the end of next 
section. 
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RANDOM LIMITATIONS 


Thus far, our analysis has assumed that both T(p^) and Its Inverse 
T(p) to be noise free* In fact, the actual measurement Is random signal 

A A 

X(Pq) which differs from T(p^) because of the contribution of noise. We 

AAA 

shall use the notation X * T + n and X ■ T + n to denote the noisy 

A 

measured and reconstructed profiles, respectively, where n and n denote 
the corresponding additive noise components. Because of the random nature 
of X, the desired profile T will never be exactly known, and It Is our 
objective In this section to define (for any given X) a confidence interval 
within which the real and imaginary components of T fall with some 
a priori specified probability, "^he tradeoff between the width of these 
confidence intevals and resolution Is examined and demonstrated with an 
example. 

Our results will be formulated in terms of the opacity t and the 
coherent phase shift defined in terms of the real, X|^, imaginary, 

Xj, and magnitude, Xj^^, components of X by the two equations 

(25) 


and 


= tan'k*|,/X,3 (26) 

Both t and relate directly to the parameters of the physical model 
assumed for the rings (Marouf et al. , 1982). We emphasize here again that 
T and as given by Eqs. (25) and (26) are only statistical estimates of 
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their exact (unkonwn) values t* and ♦ * (obtained by replacing all X's 
by T‘s in the previous two equations). Confidence intervals for x and 
are obtained from the corresponding confidence "region” for Tj^ and Tj. 


Noise Model 

A 

The dominant contribution to the additive noise component n comes from 
thermal noise generated at the front end of the Earth-based receiving station 
(Eshleman et al. » 1977). The two quadrature components of this noise may be 
adequately modeled by uncorrelated band-limited (bandwidth = B) gaussian 
random functions with zero means and identicl autocorrelation functions, with 
the total noise power equally divided among the two con^onents. 

To obtain estimates of the noise power P|yj we used the complex samples 
of the coherent signal measured in a bandwidth B = 18750 Hz during time 
periods before and after Voyager 1 was occulted by the rings. These "free 
space” samples were then used to compute estimates of the mean and variance 
over equal time intervals T^^. After normalizing the mean value to unity (as 
required by our definition of transmittance), the sum of the variances of the 
two quadrature components was taken to be the desired estimate of P^ The 
results slightly fluctuate around the estimate 



(27) 


where b = 5.02 and 3.55 for the 3,6 and 13 cm wavelength signals, 
respectively, and the second equality expresses T^ in terms of the radial 
distance Ap^ traveled in the ring plane and the radial veloctiy component 
p*. a convenient representation for later use. Clearly, we assume p* to 
be nearly constant over the period T^. 
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In all calculations to follow we assume to be a known quantity whose 

A 

magnitude Is given by the estimate Pj^ of Eq. (27), A check of the accuracy 
of such assumption may be obtained from the following formula for the 

A 

normalized mean square error In Pjj^ (Bendat and Plersol, 1971, p. 176) 


VarCP^] 



1 


(28) 


Thus, fcr example, for — 100 m and p’ = 60-80 km/s (see Table I), 

T^ i 1.25 ms, which when used In Eq. (27), together with B = 18750 Hz, 
gives the error In P|^ to be less than about 20%; for Ap^ 1 1 km, on the 
other hand, the error Is less than about 7%. 

We have also carried out numerical evaluations of the signal variance 
during occultatlon by several regions of the rings. Such variance preswnably 
include contributions from the scattered signal which, as far as the coherent 
signal 1$ concerned, constitutes an additional gausslan noise component. For 
all integration times considered ( ~ scattered signal 

A 

contribution to Pj^ was negligible and was certainly within the uncertainty 
expressed by Eq. (28). 

Note that the statistical characterization described above applies to the 

A 

noise n before being processed through the inverse Fresnel transform (6). 

The linearity of the transform preserves the gaussian nature of the noise, and 
one can show that the quadrature components of the resultant noise n remain 
uncorrelated and carry equal shares of the total noise power. The noise power 

A 

is the same as P^ in Eq. (27) provided that the window width Is chosen so 
that * Ry, i.e., W = 2F^/Ap^ (c.f., Eq. (10)). 

One noise component not included in our model is the random fluctuations 
In the phase of the reference oscillator, which Introduces multiplicative 
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noise. Although presumed to be negligible, further analysis is needed to 
confirm such assumption. 

Confidence Intervals 

With measurement noise now Included, consider the complex plane whose 
orthogonal axes are the unknown variables T^. and Tj (see Fig. 5), A given 
pair (X[^,Xj), obtained from the inverse transform, is represented by a point 
in this plane which, because of the noise, is located at a position that 
differs from the actual pair (Tr,Ti). One cannot locate (T|^, Tjj exactly, but 

can inquire in a statistical sense about the region in this plane that 
contains (Tr,Tj) with some specified probability given knowledge of (X|^,Xj). 
This is a classical problem in statistical inference (v. ,e.g. , ? ), whose 

answer is obtained by noting that the random variable C, defined by the 
equation. 


C = 2[(X^ - + (Xj . (29) 

has a X - distribution of 2 degrees of freedom. For a given (Xj^,Xj), Eq. 

(29) describes circles in the (Tj^, Tj) plane whose peripheries are contours of 
constant values of C. The particular circle that contains the pair (Tj^.Tj) 
within it with some chosen probability (1 - <*) defines the so called (1 -a)% 
confidence region. The radius of such circle is /P^C^/2 , where is the 

value of C such that Prob[C > C^] = a. Tables of the cumulative 
X^-distribution gives for the 50, 90, and 95% confidence regions to be 
1.39, 4.61, and 5.99, respectively. Clearly, the region becomes larger as 
higher confidence is desired and/or as the noise power P|»^ increases. 
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From Fig. 5, one may also directly Infer (1 - a)% confidence Intervals 
for T and using the dashed contours tangent to the circle In that 
figure. Our results will be expressed In terms a "threshold opacity" 
defined as 


■ -“0 ’ 


(30) 


which corresponds to the threshold case when the radius of the circle In Fig. 
5 Is such that the origin point is located on the circle's periphery. For a 
measured pair (x,^ ), the (1 - a)% confidence Intervals that simultaneously 
contain the actual pair (x*,^^*) can be shown to be Xj^ i. x* 1 x^^ and 

l^c - 4>c*! - ^♦c* 




'U ' - ^'*0 




(t-XTH)/2>‘c 


= s1n“ [e 




+ 13 

0 < X < • 

(31a) 

- 1] 

0 « t < Tjh 

(31b) 





0 < T < tTH 

(32) 


V 

— 1 
z 



Measured values of x close to or greater than x^j^ correspond to situations 
In which the phase has little or no meaning and only a lower bound on th.^ 
actual opacity x* can be established. 

As obvious from Eqs, (30) and (27), factors affecting x^^^ are the 
occultatlon geometry (Uq,p*), the desired resolution (^Pq). and the chosen 
confidence level (C^). For Voyager 1 geometry = 0.103, and p* in Table 
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I), F1g. 6 shows variation of with resolution Ry{Ry “ for the two 

case o « 0.5 and 0,1 (50 and 90% confidence); solid lines correspond to p^* 
in inner ring C while dashed lines correspond to outer ring A. Plots for the 
3.6 and 13 cm wavelengths are shown. Note the degredation in values of 
at the longer wavelength, mainly beause of approximately 15 dB free-space SNR 
difference between the two signals. 

Fig, 7 shows an example of the varitions in the uncertainty intervals 
with measured opacity t and data processing resolution Ry^. The example 
corresponds to the choices a * 0.5 * 1.39), p* =83 km/s (outer ring A), 

and X = 3.6 cm. Similar families of curves for other choices can be 
generated from the results above. The intersection of a vertical line drawn 
at a measured opacity x with the curves labeled by the same resolution 
number give (in cycles) from the upper set of curves and x^ and 

from the lower ones. Note, for any given resolution, the rapid divergence in 
the width of the x-confidence interval as x approaches Xj^ and the 
"tightening" of the interval and increase in as resolution is degraded. 

In the light of the above results, a comment about the resolution limits 
in Table II is in order. Because these limits were derived assuming unlimited 
SNR, they clearly apply only to regions whose opacity x is small compared to 
the corresponding x^^. If, on the other hand, x is close to or larger than 
the confidence intervals become intolerably large to allow any 
meaningful quantitative interpretation of x and (j»^. To narrow down the 

uncertainty, one is forced to perform coherent data averaging over longer 
intervals Ap^ in order to reduce P|ij, hence increase Xy^^ The price paid, 
of course, is degredation of resolution. As obvious from the 3.6 cm-X curves 
in Fig. 6, Xy^^ for the Voyager 1 geometry is less than about 0.8 for 
Ap^ i 1 km, mainly because of nearly closed rings at the time of occultation 
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(B_ ■ 5.92®). Thus, over a good fraction of the rings, the resolution will be 
noise limited. One must keep In mind, however, that this conclusion applies 
only to tne particular geometry of Voyager 1 at Saturn. To highllte this 
point, we show In Fig. 6 values of corresponding to a hypothetical 
Voyager>11ke occultatlon experiment with the rings assumed to be fully open (B 
■ 27®; the vertical scale on the right side of Fig. 6). A 100 meter 
resolution nearly everywhere within the rings (radio opacity i 1.5) Is 
certainly within the realm of such experiment. A nearly Ideal target In that 
regard Is, of course, the nearly fully open (£ “ 90°) Uranian rings, 
expected to be occulted by Voyager 1 In January 1986. 


To be Completed 
Example I Revisited 

To Illustrate: sharpness of the Inner edge of the Encke Division at 

100 m resolution; implications for edge thickness; 
detached ringlets within the gap; wave phenomena 
Interior to the Inner edge. 

Other Possible Examples to Choose From 

Example II - Ring F 

To Illustrate: Why 100m resolution is necessary for narrow features 

(bifurcation of ring F); implications for rings of 
Uranus. 
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Example HI > Mimas 5:3 Density and Bending Waves 


To illustrate: Phenomenon masked by diffraction and revealed by the 

inversion. 

Example IV - A condensation in ring C 

To illustrate: utility of recovered phase & opacity at two 
observation wavelengths. 

Example V - Maxwell's Gap (ring C) 

To illustrate: sharpness of edges of embedded ringlets; fuzziness of 

edges of the mother gap; tradeoff between resolution 
and SNR. 

Example VI - Some Gaps in Cassini Division 

To illustrate: sharp outer edges; diffused inner edges; similarities 

among two gaps; possible wave phenomena next to a 
gap. 


CONCLUSIONS 
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